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ABSTRACT 

A connection of a variety of tight-binding models of noninteracting electrons 
on a rectangular lattice in a magnetic field with theta functions is established. 
A new spectrum generating symmetry is discovered which essentialy reduces 
the problem of diagonalization of these models. Provided that one knows 
one eigenvector at one point in the parameter space of the corresponding 
Harper equation one knows an eigenfunction of the corresponding model in 
the whole range of momentum singlet out by the Landau gauge. 
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The classical problem of Bloch electrons in magnetic fields has been stud- 
ied by many authors |l| , however, despite some general rigorous results about 
the spectrum of the problem [[TJ, ^[ the analytic expressions for the energy 
spectrum and wave functions are still unknown. In what follows a whole va- 
riety of tight-binding (t — b) models of noninteracting electrons in a magnetic 
field is discussed. We shall start with a rational magnetic field, a — p/q, 
describe symmetries of these models and establish a connection with theta 
functions with chracteristics ||. In our construction the whole symmetry of 
the Hamiltonian oit — b models is used which is the invariance under contin- 
uos magnetic translations (MT's), i.e., under an action of the full Heisenberg 
group - a fact which has been ignored in previous discussions of the problem 
[|I|, The use of theta functions is then natural since they arise in connec- 
tion with representations of finite dimensional subgroups of the Heisenberg 
group H which are relevant for a description of the spectrum of the Bloch 
electron. The ordinary theta functions have appeared in the study of a re- 
lated problem || , however, so far as we know nobody applied theta functions 
with characteristics to the study of the spectrum of t — b models ||. The 
problem of diagonalization of these models or the problem of solving a poly- 
nomial equation [1]] is then reduced to one point of the Brillouin zone (BZ) 
for the momentum singlet out by the Landau gauge. This is shown by using 
modular transformations, SL(2,Z) H]. 

Let us start with some definitions. In what follows the Landau gauge 
A = B(— y,0,0) and a substrate potential periodic under translations by 
a~[ and a~2 in z and y direction, respectively, are assumed. A truncation of 
the Fourier cosine expansion of the substrate potential accomapanied with 
the Peierls substitution then leads on the Hamiltonian of a t — b model [|IJ 
including in general terms describing j-th-nearest-neighbour hopping up to 
J < n, 

H = t 1 (S a - 1 +S* i +S a - 2 + St) + ..., (1) 

where S± -j = S%) is a shift operator, S ±a - j = e ± i^ m ^, v j} j = 1,2, 

being components of the standard velocity operator, and ... in ([j]) stands for 
integer powers j of the shift opertors multiplied by a corresponding overlap 
integral tj and describing next-nearest-neighbour hopping, etc. Action of 
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S±a- :j on a function if)(r) is given as follows, 

S± a -Mr) = e ± s a i(pi+f^) ? /;(r) = ip(f± a 1 )e ±2niay / a2 , 

where according to Ref. a = $/$ D) & and $ G = hc/e being the magnetic 
flux through an elementary plaquette and the flux quantum, respectively. 
Although recently a paper appeared where the use of the Peierls substitution 
is shown to lead to incorrect conclusions for total energy of electrons in a 
magnetic field ||, single electron energies are not affected (after some shift 
and rescaling) by this substitution . The displayed part of the hamiltonian 
Ti which is the Hamiltonian of the nearest-neighbour t — b model leads on 
the so-called Harper equation (HE) [|l|], 

(Hu)(n) = u(n + 1) + u(n - 1) + 2 cos(27rcm + ®)u(n) = Eu{n). (3) 

Here u(n) = g(n), where g(n) enters a parametrization of the wave function: 
ip(ma\,na2) = e lkim g(n), k\ being component of the Bloch momentum, and 
6 = kia%. 

Commutation relations for the components of the velocity operator in a 
magnetic field imply the following commutation relations for the shift oper- 
ators 5± -j, j = 1, 2, 

Sa^Sajz = Sd^Sa{G ■ (4) 

By a simple consideration one finds that 1/a roots (powers) of the shift 
operators, i.e., the Azbel operators A and B |, i := S^ a , B := S^( a , 
commute with the Hamiltonian 7i. In contrast to the shift operators their 
commuation relation is 

AB = BAe 2ni/a . (5) 

The above defined operators are not the only ones which commute with TC. 
One can choose for generators of operators commuting with 7i components of 
the velocity operator transformed under spatial parity. Indeed working in the 
above Landau gauge one can check that any of the two operators from the set 
{pi + - c By,p2\ commutes with both operators from the set {pi, p2 + -Bx} 
and vice versa. Hence, the operators T±^, j = 1,2, called operators of 
magnetic translations (MTO), 

T, ■— p^Ti^P 1 T, .— P ±i{j-a 2 p2+2max/a 1 ) tc\ 

± ±t,a\ ■ c j ± ±k,a,2 • ° 5 V u / 
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commute with 7i for any £. One can justify that the spatial parity trans- 
formed components of the velocity operator are not gauge equivalent to the 
primary ones (the model under consideration is not (except for a = 1/2) 
parity invariant), and hence MTO's are independent of the shift operators, 
as well as of A and B. Whenever £ is an integer we shall call MTO's as 
the lattice or integer MTO's. The lattice MTO's form a projective (ray) 
representation of the translation group. For any a: 



Therefore, in virtue of flg,0), as a maximal set of commuting operators can 
be taken any set of the form 

{H, A pi ,B P2 ,T aid{ ,T a2d2 I pj E Z, p 1 p 2 ~ a; Oj e R,(J\(J2 ~ 1/a; j = 1,2}. 

(8) 

In a rational magnetic field a = p/q, p and q being relative prime inte- 
gers, one can always set p\ — 1, p 2 — p, &i — 1, <J 2 = g. Note that 
the "minimal" commuting MTO's are T a -j and T qd - 2 / p in this case. However, 
Tqa 2 /p doesn't relate points which are related by 7i and hence is unimpor- 
tant for the classification of the spectrum of 7i. As a whole t — b models 
have more symmetry than the model described by the Hamiltonian H = 
(l/2m)(p— ^A) 2 + V(r), with a general substrate potential V(x, y) periodic 
under translations by dj. TC is invariant under continuos MT's (if consid- 
ered in a suitable Hilbert space (see below)) while H only commutes with 
the lattice MTO's || which imply the following periodicity conditions for an 
eigenfunction if)(r) of H, 

fx ±a u y) = fix, y), fx, y±a 2 ) = fix, y ) e ^ iax / a \ (9) 

f and f being from the same degenerate set as ip. 

Moreover, apart continuos MTO's, H does not commute with the opera- 
tors B and A. Therefore the procedure which leads on the t — b Hamiltonian 
7i give a more symmetric Hamiltonian than the primary one. This is why, 
as we shall see, the problem of finding the spectrum of 7i for a given crystal 
momentum k and a given magnetic field can be essentially reduced. How- 
ever, a crucial difference between finite-order differential equation and finite- 
difference equation (i.e., infinite-order differential one) should be noted. The 
additional symmetries of 7i with regard to H don't lead to any additional 
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degeneracy of energy levels since they don't relate points which are related 
by the Hamiltonian. The additional symmetries result in that we have a 
continuous family of equivalent Hilbert spaces associated with a fixed lat- 
tice. They differ each other by the boundary conditions imposed (see below; 
for a similar case see 0). Therefore, in a rational magnetic field a = p/q the 
spectrum of either 7i or H can be classified by irreducible representations 
of the magnetic group or by a discrete subgroup of the Heisenberg group 
generated by {T a -,T ?( g}, respectively []T0j| . Action of T a - 2 which is not in a 



commuting set @ then leads to a g-fold degeneracy of energy levels which 



implies the well-known Diophantine equation for the Hall conductance a JTI 
pa + qm = 1, m being an integer. Irreducible representations of the mag- 
netic subgroup are one dimensional and can be classified by quasi-momenta k 
which take on the values in the first magnetic BZ (MBZ). Since the spectrum 
of either H or 7i is periodic under translations a — > a;± 1, it will be assumed 
that a G [0, 1). 

The Heisenberg group G, which is nothing but the central extension of 
the group of ordinary lattice translations, can be defined as a set of elements 
{(A, a, b)\X G C*;a,b G R} with the multiplication law: (A, a, b) (A', a', b') = 
(\\'e~ 2m ( ba '~ b ' a \ a + a', b + b'), where C* is a unit circle in the complex plane. 
Let us denote by qT a discrete subgroup of G, qT : = {(1, qa, b)\a, b G Z}. We 
shall denote by V the Hilbert space of entire functions f(z) with the norm 
induced by the scalar product, 

= J exp(—2iry 2 / IrriT)\f{x + iy)\ 2 dx dy, (10) 

where the integral is taken over the elementary periodicity domain, r being a 
modular parameter (modulus). Let V q be a subspace of V invariant under qT. 
The action of (1, qa, b) G qT on f(z) G V q is given as usual, (1, qa, b)f(z) = 

e nib*T+2nib( Z + g a) f ( z + qa + bT ^ and £ ^ if and Qnly if = J2 ne{1/q)z C n 

exp{irin 2 T + 27rmz}, with c n+ i = c n , i.e., V q is a q— dimensional (complex) 
subspace of V. Then a discrete subgroup G q , G q := {(A, a, b)\X G C*,a G 
Z,b G (l/g)Z}/(mod qT) = C* x Z/qZ x (l/q)Z/Z, C* being the cyclic 
group of q— roots of 1, commutes with qT. Following readily arguments in 
[§] step by step one can show that the finite group G q acts irreducibly on V q . 
Moreover, one has even an analoque of the Stone-von Neumann theorem for 
discrete subgroups of the Heisenberg group . Because of irreducibility the 
action of G q on V q determines a canonical basis for V q and G q acts in a fixed 



5 



way. The standart basis of V q is given in terms of theta functions with upper 
characteristic £a (modulo a constant), where £ — 0, 1, . . . , q— 1 0. For a given 
modulus t let the complex variable z be defined by z = z(t) := ray/a2 and 
let us consider xja\ as the lower characteristic of the modified theta function 

£a 

g[x/ai](r\roty / 02) defined as follows, 

g[X](r\Tay/a 2 ) := e mTy ^ la h\,%,\{r\z{r)). (11) 

0[»/oi](r|ra y) E V a is the usual Jacobi theta function with characteristics 
I, 

0[x/a 1 ]{r\ray/a2) = exp{7rir(n + £a) 2 + 27ri(n + £a)( \-ra — )}. (12) 

n=-oo G l fl 2 

In what follows we shall rescale our coordinates according to: Xj/aj — > Xj. 
In the following step we shall introduce the momenta (fci,^) and define 
functions g klk2 p](r\Tay) := e ifclX+ifc2 ^[ £ "](r|TOy). Let us denote the Hilbert 
space of functions g klk2 fx] with the scalar product ( ]T0| ) by IV a . For a given 
modulus r the lattice MTO's (considered as a subgroup of the Heisenberg 
group) provide the unitary irreducible representation (UIR) of G q in W a , 
G q 3 (1,±1,0) ^ T±a{, G q 3 (1,0, ±a) T ±a -. Note that in a rational 
magnetic field, a = p/q, irreducible action of the lattice MTO's requires that 
the Hilbert space be specified by periodic boundary conditions (PBC) on the 
square defined by translations by q lattice spacings in both directions despite 
that already and T qd - 2 commute (0). The condition that a state if)(x, y) : = 
Yltdtg klk>2 [ x](r\ray) is an eigenfunction of the rescaled Hamiltonian TC, TC — 
Tt/h at the point (x ,y ) of the plane is tantamount to the condition that 
the g-dimensional vector d = (d a , d\, . . . , d q -i) solves the HE, 

e~ iK2 d e+1 + e iK2 d e -! + 2 cos^ + 2vr£a]d £ = ed e , (13) 

at Ki — k\ + 2iray and K 2 = k 2 — 2nax , where e — E/t\. In the matrix 
notation the vector d has to be the eigenvector of the following hermitian 
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q x q matrix: 



e iK 2 




\ e~ iK2 



-iK 2 

c x 

c iK 2 







D iK 2 






-iK 2 



iK 2 (J 



(14) 



where Cg stands for 2 cos(Ki + 27i-£a). To show this one can use the behaviour 
of g[ e x](r\ray) under lattice translations, 



g[x±i](T\Tay) = e ±2mia 



gp]{T\ra{y±l)) 



^ iax g[ {e± x 1)a }(T\ray). 



(15) 



Note that for a given parameters (Kx, K 2 ) of the HE we have a freedom to 
choose (x a , y Q ) at our convenience since a change in (x a , y Q ) can be compen- 
sated by adjusting the values of the momenta (kx, k 2 ). Due to the periodicity 
properties (pT5|) of theta functions the same HE (for a given momenta (kx, /C2)) 



( p!3|) will be repeated at the points of the lattice (x Q + qm, y Q + nq) , m and n 



being integers. Since, for a given (x Q , y Q ) 



g[x +m± q }(T\Ta(y + n)) = g[x +m\(r\Ta(y + n)) 



g[x +m](r\Ta(y + n±q)) 



J ^fiq2ira Xo 



g[x +m](r\Ta(y + n)), 



(16) 



basis functions g klk2 [ e x}(r\Tay) satisfy the same boundary conditions (BC) on 
the elementary periodicity q x q square at any point of the (x Q + m,y + n) 
lattice. This corresponds to the fact that these points are connected by the 
Hamiltonian 7i and states on this lattice belong to the same Hilbert space. 
Because of (|16D any basis function carries an internal Bloch momentum k 2 = 
—2ixax y on the (x ,y ) lattice. Thus, the state ijj(x,y) := ^ dgg klk2 [ x 
}(r\ray) carries the Bloch momenta (ki,K 2 ) on the (x ,y ) lattice. The 
Hilbert space on a given (x ,y ) lattice is determined by the scalar product 



(g k ^[n9 K H%)w^ yo 



^2 g klk2 [xo+m](T\ray + n)g^ k ' 2 [xo+rn](r\Tay + n). 

0<m,n<q—l 



(17) 
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One can check that W a (x Q , y Q ) is the Hilbert space, 7i is the hermitian oper- 
ator, and that g klk2 [x]'s are orthogonal for a given (k\, k 2 ), (. [ s "](r|ra?/)|. [S 
}(T\Tay))w a (x ,y ) = Q5 se J2 q n= 1 9[^}(T\ray )g[x^}(T\Tay ) The notion of eigen- 
function on a given lattice is selfconsistent since if ip(x, y) = J2e di(Ki, K 2 )g klk 
^x]{r\ray) is the eigenfunction of 7i at a given point of the (x ,y ) lattice 
that it remains to be eigenfunction of Ti for all points of (x Q , y Q ) lattice. In- 
deed, at the point (x ,y + 1) one obtains the HE at (K\ + 2na,K 2 ). By 
using the periodicity properties (|15D of the basis functions one can pull back 
the equation to the original point (x Q , y Q ). The above statement will be then 
proved if one shows that 4>{x,y) = Y J £^~ 2max di-i(Ki,K 2 )g klk2 ^S}{T\Tay) 
solves the HE at the above parameters at the original point (x ,y ). How- 
ever, this can be checked rather straightforwardly. Similarly, at the point 
(x + l,|/o) one has to show that tp(x,y) = J2e de(Ki, K 2 )g klk2 [ e x](r\ray) 
solves the HE at (K\,K 2 — 2ira). This is equivalent to show that ip(x,y) = 
J2ee 2niia d e (K 1 ,K 2 )g klk2 fx}(T\Tay) solves the HE at the above parameters 
at the original point (x ,y ). This can be again justified rather straightfor- 
wardly. Therefore any solution df,(K\,K 2 ) of the HE defines eigenfunction 
ip{x,y) = J2ede{Ki, K 2 )g klk2 [ x]{r\ray) of TC in the whole continuous fam- 
ily of the Hilbert spaces of states on the (x ,y ) lattices if momenta (ki, k 2 ) 
and the lattice (x ,y ) are adjusted such that the parameters (Ki,K 2 ) are 
kept fixed mod 2ira. The continuous MTO's can be viewed as operators 
which map states from one Hilbert space of states associated with a given 
lattice to another. By means of them one can put the eigenf unctions of 7i 
for different momenta (ki,k 2 ) and, of course, for different (x ,y ) lattices, 
on a given lattice, i.e., to the same Hilbert space. Let us consider a state 
ip{x,y) = Hi dtg klk2 fx\{T\Tay) e W a (x ,y ), i.e., the (x a , y a ) -lattice section 
of W a , with the components dt = de(K\, K 2 ), which solves ( |13D for (Ki, K 2 ) 
at the point (x ,y ). Then ip' := T tjd -.ip e W a (r Q + ejdj), j = 1,2, satisfies a 
different BC and it is an element of W a {r + tjdj). It ceases to be an eigenvec- 
tor since its components dg have to satisfy the HE at (Ki+2iiae 2 , K 2 —2^at\), 
i.e., 

e- i{K2 - 2nQei) d e+1 + e i(K2 - 27raei) ^_i + 2 cos^ + 2irla)dt = ed e . (18) 

Thus, if we wish to make the property that some state ip is an eigenfunction 
invariant under continuous MT's, T e . -., we have to simultaneously trans- 
form either the components df(Ki,K 2 ) according to de(Ki,K 2 ) dt(K\ + 
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27rae 2 , K 2 — 27rae\), or the momenta (ki, k 2 ) by the following rule: (k\, k 2 ) 
(fci — 2irae 2 , k 2 + 2naei). First, let us consider the combined transformation 
(the shift by T e . -. and the transformation of d e (K 1 , K 2 )) by C/(ei) and £/(e 2 ) 
separately. Therefore the MTO's don't commute in general the image of 
Y\j s U(ej s ) (s being an index, one or two, which distinguishes between trans- 
lations in x and y directions) is only defined modulo an overall phase factor. 
To get rid of it we shall define the operator Q(ei, e 2 ), 

Q( ei ,e 2 ) := [e^Cn^Je-^II^.). ( 19 ) 

3 s 3 s 

where the both products are ordered in the same manner and e s = J2j 3 £j s ■ 
From the above it follows that Q(ei,e 2 ) defines a periodic function on the 
(to, n)-lattice, to and n being integers. Similarly we shall define the other 
combined mapping R(ei,e 2 ) (shift by T e a -. and the transformation of mo- 
menta (ki,k 2 )) by i?(ei,e 2 ). The operators Q(ei,e 2 ) and R(ei,e 2 ) defines 
canonical mappings between different Hilbert spaces W a (x ,y ) under which 
eigenvectors are mapped on eigenvectors. Thus we know that for any e\ 

and e 2 both Q(e u e 2 )^ {Xo , yo) := Eed i (K 1 + 2nae 2 , K 2 -2itae 1 )g k ^ +2 ™^[x+e 1 
](r\ra(y+e 2 )) and R(e u e 2 )^ {xo , yo) := ^ e d e (K u K 2 )g k ^ +2 ™^ [ a! + ei ](r|ra(y+ 
e 2 )) are eigenfunctions of 7i at (x ,y ) if ip( Xo , yo ) :— J2e dt(K 1 , K 2 )g klk2 fx 
](r\ray) does so. By using the modular transformations we can explicitely 
write the induced change of the components di down. Let us suppose that 
there exists a relation between our basis functions at different points (x, y) 
and (x,y), say 

g k ^p]{r\ray) = M es (r\x,y)g k ^fi](f\ray). (20) 

Our main idea is as follows. Let us assume that there exists a parametrization 
of the components de = F(d~e) of a vector ip(x,y) = ^2ide,g klk2 fx]{r\Tay) 
which exactly cancels the effect of Aie s (r\x, y) in the transformed basis under 
the action of 7i. Then the condition that ip is an eigenfunction on a given 
point (x Q , y Q ) lattice is tantamount to the condition that both, di and d~£, solve 
the HE. However, the equation for de is in general at different point of the BZ 
than the equation for de. The parametrization then in turn induces a mapping 
between components of eigenvectors at these points of the BZ. In what follows 
we shall look for relations between our basis functions of the above type 
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(pop . The simplest relations between the spectrum of H at different points 
of the BZ follow from the behaviour of g klk2 [ e x]'s under lattice translations 
(ID): d t (K x , K 2 ) = e 2 ™ ea d e (K u K 2 - 2na) and d t (K x , K 2 ) = e^ iax d^x(Kx + 
2ttci,K 2 ). Since g klk2 [ e £\(T\ra(y±l/a)) = e^ 2m g k ^x\{r\Tay) the action of 
Q(0, ±l/a) implies the relation di(K\, K 2 ) = dg(Ki ± 2ir, K 2 ). Now we shall 
show that the required relation ( p0|) is provided by modular transformations 
[§]. The modular transformations are transformations of the form, 

r^r=^, ,^i = -5-, (21) 
cr + a cr + a 

the parameters a,b,c,d being integers, ad — be = 1, which form the group 
of modular transformations, SL(2,Z). If the Jacobi theta functions with 
chracteristics are transformed under the modular transformation || then 
the transformation induces the following identity between the modified theta 
functions: 

g k ^[Jbay\(r\raay) = (cr + d)-^ 2 u e e~ lnaba2y2 g klk2 f±](t\tay), (22) 

where £a = £ad—cx+cd/2, x = ax—b£a+ab/2, and ui is a phase factor which 
doesn't depend on r and y. One of the key points in our construction is that 
the (finite dimensional) Hilbert space W a constructed above is not the unique 
choice. By using the periodicity properties of modified theta functions (|15|) 
one can show that the representation space of the lattice MTO's ([?]) can be 
taken to be the Hilbert space generated by the modified theta functions from 
either right or left hand side of Eq. (|22"D, as well. Of course, the parameters 
of the modular transformation ([H]) have to be such that (a — 1), ba, and 
(d— 1) are proportional to q. Then, the contribution —Ma and ab/2 to x can 
be ignored, the upper characteristic remains rational (modulo a constant) 
and proportional to a, etc. Moreover, one can show that doesn't depend 
on I in this case. One can also check that these modular transformations are 
not empty. One of them is, e.g., a = (aq 2 + 1), d — (dq 2 + 1), b = dq 2 , and 
c = a(q 2 + k), whenever a + d = had. 

Now, if one considers a vector d = J2e d((Ki, K 2 )g klk2 [x+bay](T\raay), then 
the condition that d is an eigenvector of Ti at (x Q , y Q ) is tantamount to say 
that the components of d satisfy the Harper equation at (Kx, K 2 ) 6 BZ, where 
Kx — kx + 2n-ay , and K 2 = k 2 — 27iaa(x a + bay ). Now, if one parametrizes 
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the components d^R^K?) of d according to de(Ki, K 2 ) = e lnaba2( - y ° ^ 2 dt 
then one gets the HE for dg at (Kx,K 2 + 2iraba 2 y ). The parametriza- 

1 2 2 

tion is taken to be such that it exactly cancels the effect of e~ l7T y if 
one looks at the transformed basis (|22|). Therefore dg,{K\, K 2 + 2itaba 2 y ) = 
e~ lnaba2( - y °~^ 2 d£(Ki, K 2 ). Keeping (Ki,K 2 ) fixed and changing {ki,k 2 ) and 
(x ,y ) one can distribute the eigenvector at (K 1 ,K 2 ) e BZ throughout the 
whole range of K 2 . 

By close inspection one finds that the same construction can be done for 
allt — b models discussed here (including asymmetric ones with direction de- 
pendendent overlap, 7^ tj). This is a rather straightforward consequence 
of the invariance of Ti under the Heisenberg group. More detailed consid- 
eration of the above symmetry with regard to the diagonalization of t — b 
models on a more general lattice and a consideration of the flux phases will 
be pursueded elsewhere. 
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